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Abstract 

We investigate sequential tunneling transport through a semiconductor 
double quantum dot structure by combining a simple microscopic quan- 
tum confinement model with a Mott-Hubbard type correlation model. 
We calculate nonperturbatively the evolution of the Coulomb blockade 
oscillations as a function of the interdot barrier conductance, obtaining 
good qualitative agreement with the experimental data over the whole 
tunneling regime from the weak-coupling individual dot to the strong- 
coupling coherent double-dot molecular system. 
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By tuning the tunnel barrier between the individual dots of a voltage-biased semi- 
conductor double quantum dot system, it has recently been possible^ to observe 
the formation of an artificial double-dot molecule (with each dot as its atomic con- 
stituents) in Coulomb blockade transport experiment s&H. As the interdot tunneling 
is increased, the series of linear conductance peaks of the two individual near-identical 
changes continuously to a series of split peaks which then form a well-defined 
Coulomb blockade oscillation series with twice the individual Coulomb blockade pe- 
riod. This period doubling transition in the Coulomb blockade oscillations closely 
follows the energetics of the transition of two fully isolated dots into a single compos- 
ite dot due to enhanced interdot tunneling. This transition raises important general 
questions on how the parameters which can be uniquely defined for the isolated sys- 
tem would renormalize in the transparent composite system. For example, to add an 
electron to a single dot requires energy in excess of the intradot interaction energy 
Mil. For two isolated dots in series in a double-dot system, supplying the required 
energy uu corresponds to the addition of one electron to each of the dots. With 
the increased 'transparency' of the system due to enhanced tunneling, the energy 
Mil required to add an electron changes as the doubling in the periodicity of the lin- 
ear conductance Coulomb peaks demonstrates. The pertinent theoretical question is 
how to characterize the increase of the 'transparency' of the composite system. The 
classical capacitive charging model attributes this transition to the interdot electro- 
static coupling energy Ui2- In this model each dot is considered to contain an integer 
number of electrons. The splitting of the individual Coulomb blockade peak is pro- 
portional to the interdot coupling uu- The saturation of the splitting in the strong 
coupling limit is explained by the increase of uu to the value of uu/2 with all the 
other system parameters assumed to be constantlll. The energy Ui2 is taken to arise 
from the capacitive coupling between the two dots, which is classically determined by 
the fixed geometrical arrangement of the dots. Thus, within the classical capacitive 
charging model the fixed geometrical arrangement of the two dots in the experimental 
system provides no physical reason for the increase of Ui2 necessary to account for the 
saturation of the observed peak splitting. It has also been notedS that the classical 
capacitance values needed to quantitatively explain the observed peak splitting are 
unphysically high. 

A recent fully quantum theorySli of charge fluctuations between the two dots 
formulates the double-dot problem in terms of the dimensionless interdot conductance 



2 



per tunneling channel g = Gint/{NchG^ /h) and the number of interdot tunneling 
channels Nch- In the limit of the tunneling bandwidth being much larger than the 
intradot charging energy, each dot is treated as an infinite charge reservoir for the 
other dot. A perturbation analysis in g or in (1 — g) reproduces the weak and strong 
interdot coupling limits of the peak splitting. This elegant perturbative analysis does 
not apply in the intermediate regime where Q « g « 1. The charge fluctuation 
perturbation analysis has not yet been extended to the nonlinear transport regime, 
which is one focus of our work. 

In this paper we study the transition from degenerate one dot Coulomb blockade 
oscillation to the coherent molecular double-dot oscillation by using a 2-site gen- 
eralized Mott-Hubbard model within an extremely simple physically motivated mi- 
croscopic confinement potential describing the double-dot system. A Hubbard-type 
model of linear transport through the single-particle states of quantum dot^Hl pre- 
dicts the distinct phases in the conductance pattern characterized by an increase of 
the interdot tunneling strength t. In the strong tunneling limit the Mott-Hubbard 
insulator-metal transition opens a transmission channel through an array of quantum 
dots!. This transition is equivalent to the formation of a single composite coherent 
molecular system with a concomitant doubling of Coulomb blockade periodicity. To 
characterize the dependence of the Hubbard model parameters on the value of the 
interdot conductance g, we use a simple phenomenological step-well model for the 
confinement potential profile of the double-dot system. We obtain the exact Hub- 
bard model predictions for the zero-bias voltage limit current peak splitting and the 
threshold voltage for the onset of conduction through the double-dot with 1 and 2 
excess electrons. 

The capacitive model of the experimental circuit configuration for the double-dot 
systemllHi is shown in Fig. 1. We consider a symmetrical configuration of two identical 
GaAs dots with the same electrostatic couplings to the common back gate and to the 
bias leads, i.e. Cgi = Cg2 = Cg, Ci = C2 = C, Vgi = Vg2 = Vg, and common 
self-capacitances Cq. The interdot capacitor with the capacitance Cint provides the 
electrostatic coupling between the two dots. Following the experiments, we set an 
asymmetric bias across the system, i.e. Vi = V, V2 = 0. We express the electrostatic 
part of the free energy of the system using the classical capacitance matrixll formalism. 
In the usual final step of going into the quantum-mechanical description we replace 
the classical excess charge on a dot by the charge density operator. Then the operator 
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of the electrostatic free energy of the coupled system of the double-dot and leads is 

= un{N^ + Nl) + U12N1N2 + eV{xbiNi + x^aiVa) + eVgXg{Ni + N2). (1) 

In Eq. (1) the Hubbard parameters are expressed through the elements of the ca- 
pacitance matrix, i.e. Un = e^^, U12 = e^^^, Xg = ^{2uu + Uu), a^w = 
f (2uii), Xb2 = ^{uu), where = Co + Cg + C, and 6 = C^^ - dJ . The 
bottom of the conduction band, considered to be the same in all GaAs electrodes 
in the system, is taken as the reference. Due to quantum confinement in the dot 
the continuous conduction band for an excess quasiparticle becomes a discrete series 
of single-particle energy levels Sa where a denotes the single-particle state including 
spin. We consider spin-degenerate single-particle levels in the dot. The quasiparticles 
are allowed to tunnel between the single-particle states in the two dots with the tun- 
neling amplitude ta- By including the kinetic energy in each dot and the tunneling 
energy, we write the total free energy operator as 

Fq = F''+ Y1 ^iaC\aCia-Y^{taC\aC2(i + h.C.). (2) 
i=l,2, a a 

The indices 1, 2 denote the spatial positions of dots, iVj = Y.a^^'^iaCia is the density 
operator, where c\a{.Cia) is a creation (annihilation) operator for a quasiparticle on 
the i-th dot in a state a. 

The double-dot is isolated from the leads so that it is coupled to them only electro- 
statically and through (very weak) tunneling matrix elements t^J^. The conductance 
of the dot to the leads Giead = 0.02 e^/h (<< e^jh) is kept constant throughout. 
(The lead-dot tunneling strengths t}'^ are estimated to be 3 fieV.) The tunneling 
hybridization energy between the double-dot and the bias leads is treated as a weak 
perturbation Ht'. 

H^'^ = Y.{t];^chci,2a + h.c.), (3) 

ka 

where the index k denotes a quasiparticle state in the leads. 

Using Eqs. (l)-(3) it is possible to calculate the current through the double- 
dot system under the assumptions of sequential tunneling and weak coupling to the 
outside leads, provided that the level separation in the double-dot system is much 
less than the width of the transmission resonance. We skip the calculational details 
which follow the corresponding single dot Coulomb blockade formalismih0 (suitably 
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generalized to the double-dot system). In particular, the Coulomb gap AVgap and the 
normalized peak splitting|ll, /, defined as the ratio of the additional energy needed to 
increase the number of quasiparticles by one to its maximum (saturation) value are 
respectively given by (with the total number of particles being 1 or 2) 



AVgap 



(2 - b2){uu + gi) - V&f (nn + eQ^ + [(2 - h)' - bl]t^ 
[(2 - hy - bj] 



(4) 



/ = (un + Uu/2 + 2t- ^{uu - ^12/2)2 + 2^^) /(f//2), (5) 

where u and t are those appearing in the Mott-Hubbard model defined through Eqs. 
(1) - (3), U is the intradot interaction energy of the isolated dot, bi = Xbi — Xb2 
and 62 = Xbi + Xb2- Thus a knowledge of the Mott-Hubbard parameters uu, uu, t 
etc. allow us to obtain the complete current- volt age characteristics of the double-dot 
system. 

We construct a simple microscopic quantum mechanical model to phenomenologi- 
cally describe the double-dot system depicted in Fig. 1. Our microscopic confinement 
model is shown as an inset of Fig. 2. The model uses two identical one dimensional 
infinite hard wall potential wells to describe the two isolated dots. Each dot is rep- 
resented by a two-step well (as shown in Fig. 2) with the lower step well of width a 
representing the intradot interaction energy of the individual dot and the rectangular 
barrier of potential height Vb and width d representing the (variable) tunnel barrier 
separating the two dots. When the barrier Vb is large (e.g. Vb —>■ 00) the tunnel 
conductance is vanishingly small and the two dot system is in the uncoupled 'atomic' 
limit whereas for small Vb (e.g. Vb —>■ 0) the dimensionless 'tunnel' conductance 
approaches unity and the system is in the composite 'molecular' limit. Within our 
one dimensional confinement model the barrier of height Vb and width d approxi- 
mates the intradot constriction. Such simple one dimensional potential confinement 
models have earlier been usedS to study quantum tunneling characteristics in three 
dimensional systems. Without any loss of generality we choose d/a = 0.1 throughout 
our calculation - the exact values of d and a are, of course, unknown but the only 
relevant point for our simple model to be qualitatively meaningful is for d << a. 

We evaluate the Mott-Hubbard interaction parameters Uu and by taking the 
appropriate expectation values of the screened Coulomb interaction using the poten- 
tial confinement model. (We take the screening length as 220A°.) The short-ranged 
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part of the Coulomb interaction is assumed unscreened and is approximated by a 
delta function potential. This particular way of calculating the intra- and inter- dot 
Coulomb interaction energies is obviously not unique, and one can use alternative 
definitions for uu and Ui2, which would give somewhat different numerical values. 
We feel, however, that within our simple confinement model our physically motivated 
definitions for the Mott-Hubbard parameters is reasonable. Similarly, the hopping pa- 
rameter t can be defined in several alternative ways within our microscopic model. We 
evaluate t as t = Asas/'2 where Agas is the so-called symmetric-antisymmetric energy 
gap between the two lowest single particle energy levels in our model double-well po- 
tential. Conduction in the coupled dot system occurs through a single spin-degenerate 
quasiparticle state at the ground state energy ei of the potential well. Finally, we need 
to evaluate the interdot conductance g within our microscopic model in order to make 
direct contact between the experimentalEM (and earlier theoreticalBi) results and our 
model calculations. The interdot conductance g can of course be exactly evaluated for 
our simple one dimensional rectangular barrier model of the point contact separating 
the two dots. Because of our hard wall potential confinement model, however, the ex- 
actly calculated ('hard') conductance for a tunnel barrier of height Vb and width dis a 
rather poor approximation (even on a qualitative level) for the experimental interdot 
point contact conductance. We have, therefore, also employed a 'soft' conductance 
approximation using a WKB expression, g ^ exp{— 2[^^p^(Vfe — £i)]^^^}, which we 
believe better represents (on a qualitative level) the adiabatic confinement potential 
expected in the experimental double-dot system. We find much better qualitative 
agreement between our theory and experiment using the 'soft' (as against the 'hard') 
conductance model, which is what we will mostly present in this paper. 

Within our highly simplified microsopic model for the double-dot system, the 
gate voltage induced lowering of the interdot barrier Vf, causes the crossover from 
two isolated dots (for large V;,) of size 'a' each separated by a distance 'c?' to a single 
composite coherent double-dot of size 2a -|- d. Thus a single tunable parameter, Vb, 
controls the 'transparency' of the system and causes the transition. Instead of using 
Vb as the control parameter, however, we follow the experimental procedure of using 
the interdot tunnel conductance g (determined completely by Vb in the 'soft' and 
'hard' approximations as described above) as the control parameter in depicting our 
results. As Vb is tuned all the parameters of our model (e.g. t, ei, uu, uu, g) vary 
as known functions of Vb- We fix the individual dot size a (~ 350 nm) using the 
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experiment a]0 value of the intradot interaction energy uu ^ 230 /ieV for the isolated 
dot. 

Our calculated Hubbard model parameters (e.g. ei, t, un, uu) are shown as 
functions of the corresponding interdot 'soft' conductance g in Fig. 2. Although un 
and Ui2 approach each other as g increases, the simplicity and the inadequacy of our 
microscopic model does not produce Un = Ui2 for g = 1. This is mainly due to 
the various approximations used in calculating the overlap integrals for the Coulomb 
energy. This discrepancy is the most severe quantitative limitation of our model. 

In Fig. 3 we show our calculated Coulomb Blockade oscillations for the double- 
dot system (in the linear regime) for four values of the 'soft' interdot conductance 
g = .16, .52, .8, and .99 at T = 87 mK {kbT << uu). Note that for the sequential 
tunneling situation considered here, the width of the Coulomb blockade peaks arises 
entirely from thermal broadening. Our calculated evolution of the Coulomb block- 
ade oscillations from the degenerate single-dot oscillations (at low g) through peak 
splitting (intermediate g) to the eventual period doubling (at large g) of the Coulomb 
blockade oscillations is qualitatively similar to experimental observations (cf. Fig. 5 
of Ref. 1). To further quantify our results we show in Fig. 4 our calculated normal- 
ized peak splitting / (Eq. 5) and the Coulomb gap AVgap (Eq. 4) as functions of the 
interdot tunnel conductance g, both for 'soft' and 'hard' conductance models. The 
corresponding experimental datJ for / show considerable scatter and our results (for 
the soft model) agree with experiment. The 'hard' model, however, disagrees with 
the experimental results for reasons discussed above. We point out (and this has been 
alluded to above) that the main quantitative limitation of our model seems to be a 
weaker (stronger) dependence of both / and AVgap on the tunnel conductance g for 
small (large) values of g than seen in experiments. 

Finally, in Fig. 5 we show for g = .16, . 8 and .99 our calculated nonlinear Coulomb 
blockade transport characteristics for the double-dot system by plotting the calculated 
current (in gray scales) as a function of both the source-drain voltage and the gate 
voltage. Again, our results are in good qualitative agreement with the experimental 
with the main quantitative discrepancy arising from the inadequacy in our 
tunnel conductance value g which is higher than the corresponding experimental 
result. 

In conclusion, using a simple single-parameter (H) one dimensional microscopic 
confinement model we calculate nonperturbatively the linear and nonlinear Coulomb 
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blockade characteristics of a double-dot system as a function of the interdot tunnel 
conductance. Our results are in reasonable qualitative agreement with the experi- 
mental results which is really the best we can hope for given the great simplicity of 
our model. We can obtain better quantitative agreement with experiment by using 
additional (e.g. d and a) adjustable parameters, but we feel that to be not particularly 
meaningful within our simple model. 

This work was supported by the US-ONR. The authors thank C. H. Crouch for 
helpful discussions of her two-dot transport experiments. 
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FIGURES 



FIG. 1. The equivalent circuit of the double-dot system under study. The values of 
the capacitances Ci = C2 = 38 aF. The other parameters are defined in the text. 

FIG. 2. Variation of the Hubbard model parameters un, U12, t, and e\ with 'soft' 
interdot conductance g. The conductance g is in units of 2e^ /h. Inset: The step- well model 
as defined in the text. 

FIG. 3. The conductance g^d in units of /h through the double quantum dot system 
versus the gate voltage Vg for four values of the 'soft' interdot conductance g = .16, .52, .8, 
and .99 at T = 87 mK in the linear regime {Vgd = 10 fieV). The full thermal equilibra- 
tion in the double dot system is assumed before and after each tunneling event. We use 
£2/5^1 = 2, t2/ti = 1.2 and ^'^/t}'^ = 1.2 for the single-particle spin degenerate levels. 

FIG. 4. Calculated normalized peak splitting / (Eq. 5) and the Coulomb gap AVgap 
(Eq. 4) as functions of the interdot tunnel conductance g (top and bottom panels) for 'soft' 
and 'hard' (thick lines) conductance models. The dashed lines show the boundaries for the 
scatter of the experimental data (taken from Fig. 5 of the second paper in Ref.l). The 
top thin line in the bottom panel shows the Coulomb gap for fixed values of AVgap with 
uii = 227 iieV and uu = .11 /leV. 

FIG. 5. The nonlinear current / through the double-dot system (at T = 87 mK) is 
plotted versus the values of the gate voltage Vg and the bias voltage Vsd for three values 
of the interdot 'soft' conductance g = 0.16, 0.8, 0.99. The brightest shades in the plots 
correspond to I = 59, 76, 78 pA for the graphs from top to bottom. The states with N = 
to N = 8 electrons in the double dot system contribute to / in the shown parameter space. 
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Figure 1. R. Kotlyar and S. Das Sarma, Nonlinear Transport through Coupled 
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Figure 2. R. Kotlyar and S. Das Sarma, Nonlinear Transport through Coupled 
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Figure 3. R. Kotlyar and S. Das Sarma, Nonlinear Transport through Coupled 
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Abstract 



We investigate sequential tunneling transport through a semiconductor 
double quantum dot structure by combining a simple microscopic quan- 
tum confinement model with a Mott-Hubbard type correlation model. 
We calculate nonperturbatively the evolution of the Coulomb blockade 
oscillations as a function of the interdot barrier conductance, obtaining 
good qualitative agreement with the experimental data over the whole 
tunneling regime from the weak-coupling individual dot to the strong- 
coupling coherent double-dot molecular system. 
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By tuning the tunnel barrier between the individual dots of a voltage-biased semi- 
conductor double quantum dot system, it has recently been possible^ to observe 
the formation of an artificial double-dot molecule (with each dot as its atomic con- 
stituents) in Coulomb blockade transport experiments^^"^. As the interdot tunneling 
is increased, the series of linear conductance peaks of the two individual near-identical 
dots^^"^ changes continuously to a series of split peaks which then form a well-defined 
Coulomb blockade oscillation series with twice the individual Coulomb blockade pe- 
riod. This period doubling transition in the Coulomb blockade oscillations closely 
follows the energetics of the transition of two fully isolated dots into a single compos- 
ite dot due to enhanced interdot tunneling. This transition raises important general 
questions on how the parameters which can be uniquely defined for the isolated sys- 
tem would renormalize in the transparent composite system. For example, to add an 
electron to a single dot requires energy in excess of the intradot interaction energy 
Uii. For two isolated dots in series in a double-dot system, supplying the required 
energy Un corresponds to the addition of one electron to each of the dots. With 
the increased 'transparency' of the system due to enhanced tunneling, the energy 
Uii required to add an electron changes as the doubling in the periodicity of the lin- 
ear conductance Coulomb peaks demonstrates. The pertinent theoretical question is 
how to characterize the increase of the 'transparency' of the composite system. The 
classical capacitive charging model attributes this transition to the interdot electro- 
static coupling energy Ui2. In this model each dot is considered to contain an integer 
number of electrons. The splitting of the individual Coulomb blockade peak is pro- 
portional to the interdot coupling Ui2. The saturation of the splitting in the strong 
coupling limit is explained by the increase of Ui2 to the value of Uii/2 with all the 
other system parameters assumed to be constant^. The energy Ui2 is taken to arise 
from the capacitive coupling between the two dots, which is classically determined by 
the fixed geometrical arrangement of the dots. Thus, within the classical capacitive 
charging model the fixed geometrical arrangement of the two dots in the experimental 
system provides no physical reason for the increase of Ui2 necessary to account for the 
saturation of the observed peak splitting. It has also been noted^ that the classical 
capacitance values needed to quantitatively explain the observed peak splitting are 
unphysically high. 

A recent fully quantum theory'*'^ of charge fiuctuations between the two dots 
formulates the double-dot problem in terms of the dimensionless interdot conductance 
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per tunneling channel g = Gint / {Nchs'^ / h) and the number of interdot tunneling 
channels Nch- In the limit of the tunneling bandwidth being much larger than the 
intradot charging energy, each dot is treated as an infinite charge reservoir for the 
other dot. A perturbation analysis in (/ or in (1 — (/) reproduces the weak and strong 
interdot coupling limits of the peak splitting. This elegant perturbative analysis does 
not apply in the intermediate regime where <<(/<< 1. The charge fiuctuation 
perturbation analysis has not yet been extended to the nonlinear transport regime, 
which is one focus of our work. 

In this paper we study the transition from degenerate one dot Coulomb blockade 
oscillation to the coherent molecular double-dot oscillation by using a 2-site gen- 
eralized Mott-Hubbard model within an extremely simple physically motivated mi- 
croscopic confinement potential describing the double-dot system. A Hubbard-type 
model of linear transport through the single-particle states of quantum dots^^® pre- 
dicts the distinct phases in the conductance pattern characterized by an increase of 
the interdot tunneling strength t. In the strong tunneling limit the Mott-Hubbard 
insulator-metal transition opens a transmission channel through an array of quantum 
dots^. This transition is equivalent to the formation of a single composite coherent 
molecular system with a concomitant doubling of Coulomb blockade periodicity. To 
characterize the dependence of the Hubbard model parameters on the value of the 
interdot conductance (/, we use a simple phenomenological step-well model for the 
confinement potential profile of the double-dot system. We obtain the exact Hub- 
bard model predictions for the zero-bias voltage limit current peak splitting and the 
threshold voltage for the onset of conduction through the double-dot with 1 and 2 
excess electrons. 

The capacitive model of the experimental circuit configuration for the double-dot 
system^^"^ is shown in Fig. 1. We consider a symmetrical configuration of two identical 
GaAs dots with the same electrostatic couplings to the common back gate and to the 
bias leads, i.e. Cgi = Cg2 = C^, Ci = C2 = C, Vgi = Vg2 = Vg^ and common 
self-capacitances Cq. The interdot capacitor with the capacitance Cint provides the 
electrostatic coupling between the two dots. Following the experiments, we set an 
asymmetric bias across the system, i.e. Vi = V, V2 = 0. We express the electrostatic 
part of the free energy of the system using the classical capacitance matrix^ formalism. 
In the usual final step of going into the quantum-mechanical description we replace 
the classical excess charge on a dot by the charge density operator. Then the operator 
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of the electrostatic free energy of the coupled system of the double-dot and leads is 

= un(iV2 + N^) + m^N^N^ + eV{xbiNi + Xb2N2) + eVgXg{Ni + N^). (1) 

In Eq. (1) the Hubbard parameters are expressed through the elements of the ca- 
pacitance matrix, i.e. Un = e^^, Ui2 = e^^^, Xg = ^(2uii -|- U12), x^i = 
^(2uii), Xb2 = ^{ui2), where Ce = Co + Cg + C, and 5 = Cs^ - C^t^ The 
bottom of the conduction band, considered to be the same in all GaAs electrodes 
in the system, is taken as the reference. Due to quantum confinement in the dot 
the continuous conduction band for an excess quasiparticle becomes a discrete series 
of single-particle energy levels e„ where a denotes the single-particle state including 
spin. We consider spin-degenerate single-particle levels in the dot. The quasiparticles 
are allowed to tunnel between the single-particle states in the two dots with the tun- 
neling amplitude t„. By including the kinetic energy in each dot and the tunneling 
energy, we write the total free energy operator as 

^0 = + ^^<yc\aC^a " Xl(^«'^^l«'^2« + h.c). (2) 

8 = 1,2, a a 

The indices 1, 2 denote the spatial positions of dots, Ni = ^„ c^iaQa is the density 
operator, where c^ia{cia) is a creation (annihilation) operator for a quasiparticle on 
the i-th dot in a state a. 

The double-dot is isolated from the leads so that it is coupled to them only electro- 
statically and through (very weak) tunneling matrix elements tj^'^. The conductance 
of the dot to the leads Giead = 0.02 e'^/h (<< e'^/h) is kept constant throughout. 
(The lead-dot tunneling strengths t\'^ are estimated to be 3 fieV.) The tunneling 
hybridization energy between the double-dot and the bias leads is treated as a weak 
perturbation Ht- 

H'/ = J2{tl''c\c,,2a + h.c.), (3) 

ka 

where the index k denotes a quasiparticle state in the leads. 

Using Eqs. (l)-(3) it is possible to calculate the current through the double- 
dot system under the assumptions of sequential tunneling and weak coupling to the 
outside leads, provided that the level separation in the double-dot system is much 
less than the width of the transmission resonance. We skip the calculational details 
which follow the corresponding single dot Coulomb blockade formalism®^^° (suitably 
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generalized to the double-dot system). In particular, the Coulomb gap AVgap and the 
normalized peak splitting^, /, defined as the ratio of the additional energy needed to 
increase the number of quasiparticles by one to its maximum (saturation) value are 
respectively given by (with the total number of particles being 1 or 2) 



AVgap 



(2 - b2){un + £i) - ^bljun + eiY + [(2 - 62)^ - bl]f^ 

[(2 - - hi] 



(4) 



/ = (till + tii2/2 + 2t- ^(un - ui2/2)2 + 2t2^ /([//2), (5) 

where u and t are those appearing in the Mott-Hubbard model defined through Eqs. 
(1) - (3), U is the intradot interaction energy of the isolated dot, bi = x^i — Xb2 
and 62 = + Xb2. Thus a knowledge of the Mott-Hubbard parameters Un, U12, t 
etc. allow us to obtain the complete current- volt age characteristics of the double-dot 
system. 

We construct a simple microscopic quantum mechanical model to phenomenologi- 
cally describe the double-dot system depicted in Fig. 1. Our microscopic confinement 
model is shown as an inset of Fig. 2. The model uses two identical one dimensional 
infinite hard wall potential wells to describe the two isolated dots. Fach dot is rep- 
resented by a two-step well (as shown in Fig. 2) with the lower step well of width a 
representing the intradot interaction energy of the individual dot and the rectangular 
barrier of potential height Vb and width d representing the (variable) tunnel barrier 
separating the two dots. When the barrier Vb is large (e.g. Vb — > 00) the tunnel 
conductance is vanishingly small and the two dot system is in the uncoupled 'atomic' 
limit whereas for small Vb (e.g. Vb — > 0) the dimensionless 'tunnel' conductance 
approaches unity and the system is in the composite 'molecular' limit. Within our 
one dimensional confinement model the barrier of height Vb and width d approxi- 
mates the intradot constriction. Such simple one dimensional potential confinement 
models have earlier been used^^ to study quantum tunneling characteristics in three 
dimensional systems. Without any loss of generality we choose d/a = 0.1 throughout 
our calculation - the exact values of d and a are, of course, unknown but the only 
relevant point for our simple model to be qualitatively meaningful is for d « a. 

We evaluate the Mott-Hubbard interaction parameters Un and U12 by taking the 
appropriate expectation values of the screened Coulomb interaction using the poten- 
tial confinement model. (We take the screening length as 220 A°.) The short-ranged 
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part of the Coulomb interaction is assumed unscreened and is approximated by a 
delta function potential. This particular way of calculating the intra- and inter- dot 
Coulomb interaction energies is obviously not unique, and one can use alternative 
definitions for Un and U12, which would give somewhat different numerical values. 
We feel, however, that within our simple confinement model our physically motivated 
definitions for the Mott-Hubbard parameters is reasonable. Similarly, the hopping pa- 
rameter t can be defined in several alternative ways within our microscopic model. We 
evaluate t as t = where Asas is the so-called symmetric-antisymmetric energy 

gap between the two lowest single particle energy levels in our model double-well po- 
tential. Conduction in the coupled dot system occurs through a single spin-degenerate 
quasiparticle state at the ground state energy ei of the potential well. Finally, we need 
to evaluate the interdot conductance g within our microscopic model in order to make 
direct contact between the experimentaP^"^ (and earlier theoretical'*'^) results and our 
model calculations. The interdot conductance g can of course be exactly evaluated for 
our simple one dimensional rectangular barrier model of the point contact separating 
the two dots. Because of our hard wall potential confinement model, however, the ex- 
actly calculated ('hard') conductance for a tunnel barrier of height Vb and width 0? is a 
rather poor approximation (even on a qualitative level) for the experimental interdot 
point contact conductance. We have, therefore, also employed a 'soft' conductance 
approximation using a WKB expression, g exp{— 2[ ^™2^ (V;, — £1)]*''^}, which we 
believe better represents (on a qualitative level) the adiabatic confinement potential 
expected in the experimental double-dot system. We find much better qualitative 
agreement between our theory and experiment using the 'soft' (as against the 'hard') 
conductance model, which is what we will mostly present in this paper. 

Within our highly simplified microsopic model for the double-dot system, the 
gate voltage induced lowering of the interdot barrier Vb causes the crossover from 
two isolated dots (for large V;,) of size 'a' each separated by a distance '0?' to a single 
composite coherent double-dot of size 2a -\- d. Thus a single tunable parameter, Vb, 
controls the 'transparency' of the system and causes the transition. Instead of using 
Vb as the control parameter, however, we follow the experimental procedure of using 
the interdot tunnel conductance g (determined completely by Vb in the 'soft' and 
'hard' approximations as described above) as the control parameter in depicting our 
results. As Vb is tuned all the parameters of our model (e.g. t, ei, Un, U12, g) vary 
as known functions of Vb- We fix the individual dot 350 nm) using the 
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experimental^ value of the intradot interaction 230 /ieV for the isolated 

dot. 

Our calculated Hubbard model parameters (e.g. ei, t, Un, U12) are shown as 
functions of the corresponding interdot 'soft' conductance g in Fig. 2. Although Un 
and U12 approach each other as g increases, the simplicity and the inadequacy of our 
microscopic model does not produce Un = U12 for g = 1. This is mainly due to 
the various approximations used in calculating the overlap integrals for the Coulomb 
energy. This discrepancy is the most severe quantitative limitation of our model. 

In Fig. 3 we show our calculated Coulomb Blockade oscillations for the double- 
dot system (in the linear regime) for four values of the 'soft' interdot conductance 
g = .16, .52, .8, and .99 at T = 87 mK {kbT << Un). Note that for the sequential 
tunneling situation considered here, the width of the Coulomb blockade peaks arises 
entirely from thermal broadening. Our calculated evolution of the Coulomb block- 
ade oscillations from the degenerate single-dot oscillations (at low g) through peak 
splitting (intermediate g) to the eventual period doubling (at large g) of the Coulomb 
blockade oscillations is qualitatively similar to experimental observations (cf. Fig. 5 
of Ref. 1). To further quantify our results we show in Fig. 4 our calculated normal- 
ized peak splitting / (Fq. 5) and the Coulomb gap AVgap (Fq. 4) as functions of the 
interdot tunnel conductance (/, both for 'soft' and 'hard' conductance models. The 
corresponding experimental data^ for / show considerable scatter and our results (for 
the soft model) agree with experiment. The 'hard' model, however, disagrees with 
the experimental results for reasons discussed above. We point out (and this has been 
alluded to above) that the main quantitative limitation of our model seems to be a 
weaker (stronger) dependence of both / and AVgap on the tunnel conductance g for 
small (large) values of g than seen in experiments. 

Finally, in Fig. 5 we show for g = .16, . 8 and .99 our calculated nonlinear Coulomb 
blockade transport characteristics for the double-dot system by plotting the calculated 
current (in gray scales) as a function of both the source-drain voltage and the gate 
voltage. Again, our results are in good qualitative agreement with the experimental 
data^ with the main quantitative discrepancy arising from the inadequacy in our 
tunnel conductance value g which is higher than the corresponding experimental 
result. 

In conclusion, using a simple single-parameter (V;,) one dimensional microscopic 
confinement model we calculate nonperturbatively the linear and nonlinear Coulomb 
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blockade characteristics of a double-dot system as a function of the interdot tunnel 
conductance. Our results are in reasonable qualitative agreement with the experi- 
mental results which is really the best we can hope for given the great simplicity of 
our model. We can obtain better quantitative agreement with experiment by using 
additional (e.g. d and a) adjustable parameters, but we feel that to be not particularly 
meaningful within our simple model. 

This work was supported by the US-ONR. The authors thank C. H. Crouch for 
helpful discussions of her two-dot transport experiments. 
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FIGURES 



FIG. 1. The equivalent circuit of the double-dot system under study. The values of the 
capacitances Ci = C2 = 38 aF. The other parameters are defined in the text. 

FIG. 2. Variation of the Hubbard model parameters un, U12, t, and £1 with 'soft' 
interdot conductance g. The conductance g is in units of 2e^//i. Inset: The step-well model 
as defined in the text. 

FIG. 3. The conductance g^d in units of e'^/h through the double quantum dot system 
versus the gate voltage Vg for four values of the 'soft' interdot conductance g = .16, .52, .8, 
and .99 at T = 87 mK in the linear regime {Vgd = 10 //eV). The full thermal equilibra- 
tion in the double dot system is assumed before and after each tunneling event. We use 
£2/51 = 2, ^2/^1 = 1-2 and t^'^ /t\'^ = 1.2 for the single-particle spin degenerate levels. 

FIG. 4. Calculated normalized peak splitting / (Eq. 5) and the Coulomb gap AVgap 
(Eq. 4) as functions of the interdot tunnel conductance g (top and bottom panels) for 'soft' 
and 'hard' (thick lines) conductance models. The dashed lines show the boundaries for the 
scatter of the experimental data (taken from Fig. 5 of the second paper in Ref.l). The 
top thin line in the bottom panel shows the Coulomb gap for fixed values of AVgap with 
Mil = 227 i^eV and M12 = .11 i^eV . 

FIG. 5. The nonlinear current / through the double-dot system (at T = 87 niK) is 
plotted versus the values of the gate voltage Vg and the bias voltage Vsd for three values 
of the interdot 'soft' conductance g = 0.16, 0.8, 0.99. The brightest shades in the plots 
correspond to / = 59, 76, 78 pA for the graphs from top to bottom. The states with = 
to = 8 electrons in the double dot system contribute to / in the shown parameter space. 
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